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The dwell time for dissipative quantum system is shown to increase with barrier width. It clearly pre- 
cludes Hartman effect for dissipative systems. Here calculation has been done for inverted parabolic 
potential barrier. 
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I. INTRODUCTION 



In the recent years Hartman's work on tunneling time 
has drawn much attention of the community At the 
time of it's publication, little attention was paid to Hart- 
mans theoretical work on tunneling time of wavepack- 
ets in the sixties. Hartman analyzed the temporal as- 
pects of tunneling by writing down solutions of the time- 
dependent Schrodinger equation in terms of a superpo- 
sition integral over stationary state solutions weighted 
by a Gaussian momentum distribution function. With- 
out explicitly evaluating the integrals he could infer cer- 
tain properties of the transmitted wave packet by ex- 
amining the magnitude and phase of the integrand. His 
main striking result was that under certain circumstances 
(opaque barrier) the tunneling time is independent of 
barrier length and the traversal time can be less than the 
time that would be required to travel a distance equal to 
the barrier length in vacuum. Many physicists hesitated 
to deal with Hartmans results since a very fast tunneling, 
or a zero tunneling time holds a serious consequence of vi- 
olating Einstein's postulate of Special Theory of Relativ- 
ity. Hartman effect was first re-examined by Fletcher [2| 
within stationary phase method for quasi-monochromatic 
non-relativistic particles tunneling through the potential 
barriers. Furthermore, on recent times, experiments with 
photonic band gap structures 0, 0| showed apparent su- 
perluminality. These observations as well as the theo- 
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retical predictions lead towards the phenomena of super- 
luminal barrier traversal. Regarding the explanation of 
this apparent superluminality, it is important to mention 
some publications over the past two decades or so (B-tl^j . 
Some suggestions have been made to explain faster 

than light phenomena by the concept of energy storage 
and release in the barrier region. The argument is that 
the group delay, which is directly related to the dwell 
time by an additive self interference term [l2| , is actually 
the lifetime of stored energy (or stored particles) leaking 
through both ends of the barrier. The relation between 
group delay (r G ), dwell time (t d ) and self-interference 
delay (t 1 ) is given by 

T G =T D +T I (1.1) 

When the reflectivity is high, the incident pulse spends 
much of it's time dwelling in front of the barrier as it in- 
terferes with itself during tunneling process. This excess 
dwelling is interpreted as the self-interference delay. This 
term can be successfully disentangled from the dwell time 
[l2l ] . Now if the surrounding of the barrier is dispersion- 
less, then the self-interference term vanishes .resulting in 
the equality of group delay and dwell time [l(| ■ In that 
case the dwell time will give the lifetime of energy stor- 
age in the barrier region. In a relatively recent work [l3j , 
Jakiel et.al. have shown that Hartman effect is valid 
for all known expressions of the mean tunneling time, 
in various non-relativistic approaches, with finite width 
barriers without absorbtion or dissipation. In a recent 
paper [l4| we have formulated the expression of dwell 
time in presence of dissipation by the formalism of weak 
measurement. We have shown that inclusion of dissipa- 
tive interaction precludes the zero time tunneling. But 
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in that work we have not discussed the behavior of dwell 
time with increasing barrier thickness. In this paper, it 
is our aim to calculate the dwell time theoretically for 
tunneling through a dissipative inverted parabolic bar- 
rier to find whether "Hartman effect" exists in presence 
of dissipative interaction with the environmental bath 
modes; i.e. whether the dwell time saturates with in- 
creasing barrier thickness or not. The inverted harmonic 
potential or parabolic barrier, which provides an infinite 
potential barrier can be used as a toy model for tunnel- 
ing, where exact Gaussian wave packets may be found 
as solutions. The inverted harmonic oscillator problem 
attracts a great deal of attention not only for being one 
of the exactly solvable potential in quantum mechanics 
flU - E^ but also because of having wide range of appli- 
cability in many branches of physics. It receives quite a 
formidable number of applications in many branches of 
physics from high energy to solid state theory. To incor- 
porate dissipation, we will use the formalism of Yu [23j 
and Yu et.al. [24|, where evolution of wave function and 
tunneling in dissipative system has been discussed. In 
Section II we will briefly discuss the background formal- 
ism, in Section III we will derive the dwell time using 
the discussed formalism and finally in Section IV we will 
summarize with some concluding remarks. 



II. BACKGROUND FORMALISM 

Yu et.al. have shown that for the simplest example 
of dissipative system, a harmonic oscillator coupled to a 
heat bath in a special case of ohmic dissipation, a self 
contained treatment can be used to get wave function 
evolution, where path integral technology need not be 
used. Later in another paper [23] , Yu used this particular 
formalism in the quantum tunneling case for a inverted 
harmonic oscillator to obtain the tunneling probability 
and current density. Here we will use this formalism to 
get the dwell time for such systems. Using the notations 
of Ref. [23|, [24| and quoting formulae from there, we first 
identify the total Hamiltonian of the system of inverted 
harmonic potential with a coupled infinite harmonic os- 
cillator heat bath as 



H = 
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±M(w 2 + Aw 2 )q 2 +q^2 i C j x j 
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(2.1) 



From this Hamiltonian, the quantum Langevin equation 
can be derived as 



q + r/q- u%q = f{t) 



(2.2) 



where rj is the damping constant and f(t) is the Brownian 
motion driving force given by 
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f{t) = -2_sjf I %0 COS UJjt + x j0 — — 



(2.3) 



This equation (12.21) can be solved to give 

q(t) = ai(t)qo + a 2 (t)q + [bjix j0 + b j2 x j0 ] (2.4) 

j 

where qo, qo, Xjo, Xjo are the initial position and velocity 
of the system and the bath respectively, and 
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a 2 = e 



(coshwi+ ^jsinhcjt) , , . 



6 ii = -mIo a 2 (t')cosoj(t-t')dt', 

b j2 = -ft f* a 2 {t r ) smoj(t - t')dt' ^ ' ^ 

with lo = (wg +?7 2 /4) 2 . The solution can be used to show 
that the wavefunction of the system plus bath can be 
written as 1241 



^(9, {&}. *) = iH x,-&.*) (2.7) 



Now let us assume that the initial wavefunction is a Gaus- 
sian wave packet centered at the right of the peak of the 
potential by zo has the form 

Ml-, t = 0) = (2Tra 2 y 1 ^e-^- Z0 ^/ 4a2+tkq (2.8) 
Following [H[ , the Green's function can be derived as 

,1/2 

G(q,q ;t,0)- V 2 "^ 
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Then the wavefunction ip(q,t) can be calculated as 
(2TT<j 2 )- 1 / 4 {a 1 +iu J oa 2 r 2 ) 

X g-(<?-Zo) 2 / Au 2 +i{c 2 q 2 +ciq+c a ) 



(2.9) 



(2.10) 



This is a Gaussian distribution with a width of 

cr 2 . = er 2 (a 2 + r 4 w 2 a 2 ) (2.11) 

where r = a /a and q c = aiz Q +a 2 hk/M. The coefficients 
in the exponential term are 
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C2 = ^r^ (lnae) 
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q c za 



(2.12) 
(2.13) 



dt 



e* qc^Qnai) - 2q c + ^i^a.r 2 (2.14) 
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The solution equation (I2.10j) of the time dependent 
Schrodinger equation should satisfy the expression for the 
current density 



2Me^ 



dq ) Me" 



\^\ 2 {2c 2 q + Cl ) 
(2.15) 
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Putting the value oi ip,ci,C2 we get 



d_ 

'dt [ 



(a 2 + r 4 uj 2 ai) (g-^o) 2 
1 e 2<r ,> 



J | ■_>(,/ -g c )^(ln^) + 44 c Jx^==p^e-^ 

(2.16) 

We will utilize equation (|2 . 16f) to get the dwell time in 
the next section. 



III. DWELL TIME IN DISSIPATIVE MEDIUM 
AND HARTMAN EFFECT 



discussed that the initial wavefunction is centered at the 
right of the peak of the potential, it is travelling from 
right to left. So let us now set go = 2zo and q = 0. Then 
using equation (|2.16|) we get 

J(q ) - J(g) = 4*|0nog) ( 1 + J^ 2) e-^ (3.6) 

Considering the spread of the wavefunction to be time 
dependent and replacing a be <jq in equation (|3.6p . then 
putting in equation (|3.5p . we get 



Dwell time is defined as the average number of particles 
within the barrier region divided by the average number 
entering (or leaving) the barrier per unit time. It cor- 
responds to the average time spent by a particle within 
the barrier irrespectively of whether it is finally reflected 
or transmitted. If we follow Winful's explanation 
when the surroundings of the barrier is dispersionless, it 
represents the lifetime of energy storage in the barrier 
region. The dwell time in a neighborhood of q is defined 
as the ratio between the particle number in the interval 
[q, q + dq] and the incoming current 
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where x 2 — 



a\ + r 4 uj 2 a%. Putting -^2- = (: 
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exp(-£) 



d(x 2 ) 



Again we substitute y — (,/x. Now as t 
1 and y — > £; as t — > co, x — > oo and 
equation (|3.8p becomes 



(3.7) 



(3.8) 

0,x -> 
-> 0. So 



Obviously, Eq. p.ip describes a balance equation: in 
the stationary case the injected current equals the decay 
rate of the probability in [q, q + dq]. The dwell time t d 
of a finite region within the context of a stationary state 
scattering problem is obtained via a spatial integration 
of Eq. (|3.ip . So the dwell time t d is given by 

T ° = -^ l^)| 2 ^ = Tr / P(9)dq (3.2) 

where Ji n = jj is the incident flux. We are considering 
constant incident flux for simplicity. Now we are dealing 
with non-stationary quantum states interacting with the 
environment. So the wave function and consequently the 
probability density function will also be dependent on 
time. Differentiating with respect to time: 



• d M f« dp(q,t) J 

T =Tkj qo — d(1 



(3.3) 



Using the continuity equation: 



. D M { q dJ{q,t) , M . . , . , A . 

T =-jrkl o ^ dq= Hk [JM - J{q)] (3 - 4) 

Now integrating equation (13.41) with respect to time: 



M , 

^ / [J(go) - J(l)]dt (3.5) 
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(3.9) 



The remaing problem is that the integration do not con- 
verge at y = 0. So instead of taking the limit of a; to 
be oo, we take a finitely long time t = Ti ong , for which 
x —¥ C and y — > 1. So now the dwell time is found to be 
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\- e -f+Erf{l)-Erf{0 



(3.10) 



Here l Erf stands for error function. Generally, the 
width of the potential (zq) is greater than the width of 
the wavefunction (a). Also we will study the behavior of 
dwell time with increasing width of the barrier. There- 
fore the exponentially decaying term e - ^ ' £ would be less 
significant. So neglecting this term from equation (|3.10p 
we get 



n 



8M( 



MErf(l) - Erf(0) 



(3.11) 



Expressing the dwell time in terms of potential width 
w = 2zq, we get 



2M 



D 



hk^/TT 

2M p(w\ 



V^(Erf(l)-Erf(w/2a))] 

(3.12) 



where measurements are made at t — and t = T . For 
long time measurement, we can set T — oo. As we have 



The figure shows that in presence of dissipative inter- 
action, the dwell time is not saturating with increasing 
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FIG. 1: F(f ) with the scaled potential width f 

barrier width, but it increases as the barrier width in- 
crease. Now let us compare this result with the classical 
case of a particle traveling in presence of a frictional force 
proportional to velocity (v) having the equation of mo- 
tion 

v + ^v = (3.13) 

The solution of this equation will be of the form 

v = woe" 7 * (3.14) 

Consider that the classical particle takes t cl time to 
travel the distance w c i in presence of the usual veloc- 
ity dependent frictional force. Then integrating equation 
pmp . we get 




(3.15) 



A condition must be imposed that 7 < , where tq = 
is the time taken by the particle to cross the distance 
w c i with constant velocity in absence of dissipation. 
Because otherwise the particle will lose all of it's kinetic 
energy before crossing the barrier. Now considering 7 <C 
~ and expanding the logarithmic term in equation (|3.15[) 
and neglecting 3rd and higher order terms, we get 

t cl =aw 2 cl + pwa (3.16) 

where a — 'j/vq and /? = 1/vq- From FIG.l and FIG. 2 we 
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FIG. 2: t cl with w ct taking a = .01 and j3 = .1 

see that the increment of dwell time with barrier width 
follow almost the same nature of that of a classical par- 
ticle in presence of velocity dependent frictional force. 
So here we find that inclusion of dissipative interaction 
makes the behavior of dwell time quasi-classical. Also 
such dependence of dwell time on the barrier width pre- 
cludes the afore mentioned "Hartman effect" . 



IV. CONCLUSION 

It is evident from the above analysis that due to inter- 
action with the environmental bath modes the tunneling 
time does not saturate with increment of barrier width. 
Unlike Hartman effect, in this case the dwell time de- 
pends on barrier thickness. As we have mentioned ear- 
lier, explanation of the saturation of tunneling time can 
be given by the co ncep t of energy storage and release in 
the barrier region Group delay, which is equal to 

dwell time in absence of self interference, is proportional 
to the stored energy and it saturates (in absence of dis- 
sipation), as the stored energy saturates. In presence of 
dissipative interaction, tunneling entity loses energy to 
the interacting bath modes and in the process the sat- 
uration of energy is prevented. As the barrier length 
increases, the interaction with the bath modes also in- 
crease, resulting more and more energy loss. So the dwell 
time increases with increase of barrier length instead of 
getting saturated. The continuous interaction with the 
environmental bath modes makes the behavior of dwell 
time quasi-classical. 
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